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Ph. 

' Abstract 

G^ ■ 

I We determine the masses of the light and the strange quarks in the MS'-scheme using 

our high-statistics lattice simulation of QCD with dynamical Wilson fermions. For the light 
^-H I quark mass we find m}^^{2 GeV) — 2.7(2) MeV, which is lower than in quenched simulations. 

' For the strange quark, in a sea of two dynamical light quarks, we obtain m^^"^°(2 GeV) = 

^ '. 140(20) MeV. 

o ■ 

^ ■ 

. 1 Introduction 
!>■ ■ 
G^ ■ 

"p^ ' The masses of u, d, and s quarks constitute fundamental parameters of the Standard Model. 

■ Phenomenologically, however, they remain among the most poorly known quantities within its 

' ! scenario. 

Qj I While lowest order chiral perturbation theory offers a fairly easy access to the determination 

• of ratios of quark masses from the empirical mesonic spectrumjQ], one has to apply much more 

, meticulous techniques, such as QCD sum rules[|| or lattice QCD[^, ^, in order to arrive at 

absolute values. At this stage, however, these two methods appear to lead to contradictory 
results. 

In practice, both of these approaches carry their specific merits and shortcomings. While 
the sum rule results are sensitive to the choice of parametrizations, as elaborated in ref.||5|, the 
lattice results have been established so far for pure gluon dynamics only[^, ^. 

It is therefore of considerable interest to study the dynamical effects of vacuum fluctuations, 
originating from light quarks, onto light quark properties such as their masses. This holds in 
particular for Wilson-like discretizations of fermions which, unlike staggered fermions, are free 
of flavour symmetry violations on the lattice. Until recently, however, computing resources and 
techniques were too limited to allow for the generation of reliable samples of vacuum configura- 
tions with appropriate statistics. Nevertheless, a recent rough analysis of world data|^ seems to 
suggest that unquenching from A'^^ = to Nf = 2 might have a sizeable impact on the value 
of the light and strange quark masses, lowering them by as much as 50 %. 

In this letter, we present a lattice analysis for the light and strange quark masses based on 
our measurements of the pseudoscalar and vector mesons determined in a sea of two degenerate 
dynamical quarks. We have generated three sets of gauge configurations at three sea-quark 
masses but at the same coupling; each set comprises 200 independent gauge configurations. 
Good signals in the autocorrelation functions and the use of a blocking method give us confidence 
as to the reliability of the quoted errors. 



^ including a preliminary data set from SESAM. 
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Pdyn = 5.6, Nf = 2, 16^ X 32 

Ksea 0-156 0.1570 0.1575 

number of configurations 200 200 200 

Ky — Ky Combinations 15 15 15 

/^quenched = 6.0, 16^ X 32 

number of configurations: 200 
— combinations: 15 



Table 1: Simulation details. 



The masses of the light and the strange quark are extracted from the meson data by fixing 
two mass ratios, ^ for the light and one of or ^^^^ for the strange; the lattice spacing 

is obtained from Mp. 

We identify the dynamical quarks with the degenerate doublet of isopsin symmetric quarks, 
u and d, called light quark in the following. Thus, using our data for the pseudoscalar and the 
vector at the three sea-quarks, we can extract the mass of the light quark in a sea of light quarks. 

To simulate the strange quark we need to introduce valence quarks that are unequal to the 
dynamical quarks: at each sea-quark we evaluate meson masses with strange valence quarks and 
perform an extrapolation of these masses to the physical sea of light quarks. This procedure 
allows us to calculate the masses of the K, K* and the (j) with a dynamical light quark(for K 
and K*) and strange quarks in a sea of light quarks. 

The definition of a strange quark mass in a sea of light quarks requires an analysis of lattice 
data in terms of sea and valence quark masses, for which we describe a suitable parametrization 
in section ^. As an aside, we also comment on a possible flaw in the extraction of quark masses 
from quenched data with Wilson- like fermions. 

2 Simulation details 

In table |l| we show the parameters of our simulation. In addition to the dynamical simulation we 
performed a quenched study at the matching quenched coupling^ to enable a direct comparison 
of full and quenched QCD. 

At each of our three sea-quark values, characterized by the hopping parameters Ksea; we have 
investigated zero-momentum two-point functions, 

CAB(t) = E(0|xi(0)XB(x)|0), (1) 

for the pseudoscalar and vector particles, Xps{x) = q'{x)j^q{x) and Xyi^) ~ q'{^)i^q{^)- We 
combined light-quark propagators with hopping parameters equal and different to that of the 
corresponding sea-quark, allowing for fifteen mass estimates per sea-quark. We use the gauge- 
invariant Wuppertal-smearing procedure!^ to calculate smeared-local and smeared-smeared cor- 
relators (with smearing parameter a = 4 and with 50 iterations). Both types of smearing are 
used to obtain mass-estimates by performing a simultaneous single-exponential^ fit to the data 
on time-slices 10-15. Details will be given in 0|. 

In ref.||8| we will present a detailed auto-correlation analysis. We found integrated auto- 
correlation times. Tint! for the masses to lie around 25 HMC time- units with a slight increase 

^This is the quenched /3 which yields a similar lattice spacing. 

^We have checked that two-exponential fits yield stable ground state masses. 
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towards lighter sea-quarks. We have therefore chosen to calculate propagators on configurations 
separated by 25 HMC trajectories. To determine by how many units of Tint the measurements 
need to be separated to achieve complete decorrelation we performed a blocking analysis and 
plotted the error as a function of the blocking size. At block size 6 (for /^-sca — 

0.156 and 0.157) 

and 7 (at Kgea = 0.1575) we find the jackknife errors of the error to run into plateaus^. Conse- 
quently, these are the errors we will quote in the following. A similar analysis for our quenched 
data shows no increase in error with the blocksize (quenched configurations are generated with 
an overrelaxed Cabbibo-Marinari heatbath update and are separated by 250 sweeps). 

Errors (on the blocked data) are obtained using the bootstrap procedure. They correspond 
to 68 % confidence limits of the distribution obtained from 250 bootstrap samples. 



3 Results - light and strange quarks 

The light quark mass is extracted by matching the ratioQ: 

^ = ^ = 0.1785, (2) 

using data with (degenerate) valence quarks equal to the sea-quarks. Generically, we call this 
data russ- Owing to the fact that we have but three different sea-quarks, 

"^^^^ = l(—-T^)^ (3) 



2 \ Ksea l^'< 

we use linear extrapolations in the lattice quark mass: 



"^ps,ss = c( — --^) , (4) 



^sca 

my ss = rn^^^^ + bm 



''""PS,ss 



These fits, which we call "symmetric" , are shown in figure |T[ We find the pseudoscalar mass to 
be extremely well matched by the linear ansatz (x^/d.o.f = 0.002), whereas the vector masses 
may exhibit some downward curvature (x^/d.o.f = 1.1). We find: 

= 0.15846(5) = 0.15841(5) , (5) 

= \{^-Tt)= 0-00088(6) . (6) 



m 



The corresponding lattice spacings from the rho are: 

a^^ = 2.35(6) GeV at (7) 
= 2.33(6) GeV at ^fi^*. (8) 



We use the latter value to convert to physical units in the MS'-scheme according to: 

"^ms(m) = ZM{^la)rrf^''a-^ , (9) 

with ZM{fJ-a) calculated in boosted 1-loop perturbation theory ||l^, |l3| and run the values to 2 
GeV. Throughout, we allow for a 3 % uncertainty in Zm- As a result we find: 

miig(2GeV) = 2.7(2) MeV . (10) 



*We consider this to be a conservative estimate. 

^We use the convention that massesj^ in the continuum are denoted with capital "A/", whereas lattice masses 
are written as "m" . 
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As we have outlined in the introduction, the treatment of the strange quark within the context 
of an Nf = 2 simulation requires the computation of mesons with valence quarks unequal to the 
dynamical light sca-quarks. In addition to rriss we introduce the following generic notation: 

• - one valence quark is identical to the sea-quark. 

• rrivv - neither valence quark is identical to the sea-quark. 

Furthermore, we define an effective k through ^ = | -|- -i^^, where and refer to 
valence quarks in a meson. 

In principle, we can fit mpg and mpg to independent linear functions in and ^ , 

"ips.vv = C1 + C2— + C3— H , (11) 



,2 



1,1 



"^ps,sv = c'l + 4 — + 4 — + ■ ■ ■ ■ (12) 

ftV "-sea 

However, from the requirement that all parametrizations must converge smoothly into each other 

on the symmetric line, Ksca = ny, the number of independent parameters can be substantially 
reduced. For mpc.,,,, for example, one finds C2 + C3 = c and ci = — In particular, at the 

^sca 

critical point mpg = mpg = mpg ^ = 0, valence and sea quark masses must be identical: 
'^v('^sea) = '^sea- ^his simplifies the mass equations to: 

"ipS,vv = -C-^ + C3— -F — (C-C3) , (13) 
'^sea ^sea. 

"ips.sv = -c^ + 4-^ + -;J-(c-c'3) • (14) 

"'sea "'sea "'V 
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Figure 2: Simultaneous fit of all pseudoscalar data to eq. 16. The vector data are fit to a semi- 



constrained form (see text). Symbols: * 
lattice units. 



O 



□ 



All plots in 



Defining a bare valence quark mass with respect to k^^^^ as: 
we can summarize the combined parametrization as follows ||To|| : 

2 

V "^PS,vv / 

where the variables C34 = 2c3 and C13 = 2c'^ — c have been introduced. Similarily, the vector-data 
can be parametrized as: 

rnv,sv 

\ "iV,vv 

Once the fit parameters c, C13, C34 and k^^^^ for the pseudoscalar and additionally m'^"^ and b 
for the vector particle have been determined, we can calculate all quark masses; no additonal 
fitting is needed. 

Figure |^ (left plot) shows such a combined linear extrapolation of all the pseudoscalar data 
with the ansatz of eq. The data are nicely fitted by this parametrization (x^/d.o.f = 4.4/23). 
A constrained fit to the vector-data with eq. ^ turns out to be more difficult; we find a x^/d.o.f = 
50/25, fit-parameters are given in table ^. A much better x^/d.o.f for the vector particle can 
be achieved by a semi-constrained fit in which only K^ga) ^ given in eq. is held fixed, but the 
parameters c, C13 and C34 are allowed to vary. This is shown in figure ^ on the r.h.s.. The data 
for m-sv and mw are best fit for those data points of russ which match eq. ^ well. 



2c 

C + C13 C-C13 II '^v 1 ' (16) 
C34 2c - C34 



m 



crit 



+ b 



"IpS.ss 
2 

"^PS,sv 



m 



PS,vv 



(17) 
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'^sca 


light 
"'sea 


c 


Cl3 


C34 


b 


^crit 


xVd.o.f 


symmetric fit to m-pg and niy^ss 


0.15846+ 5 


0.15841 + I 


1.99+ I 








0.328 +^g 


0.0024, 1.1 


constrained fit to eq. 16 and eq. 17 


0.15845 + 5 


0.15841 + I 


1.98+ I 


0.88 + g 


1.7+ ; 


1.02+ 5 


0.329 + g 


4.4/23, 50/25 


semi-constrained fit to eq. 17 


"'sea 


light 
"'sea 


^rho 


rho 


rho 


6 


^crit 


xVd.o.f 


0.15845 + 5 


0.15841 + I 


2.08+}° 


1.06 


2.1 




0.326 + 7 


4.4/23, 6.3/23 



Table 2: Fit results from constrained (equations 16 and 17) as well as from symmetric fits to 
the pseudoscalar and vector and the values of Kgg^ and 



Before we proceed to extract the strange quark mass, let us check for the consistency of the 
method by using the combined equation to determine the light quark mass and by comparing 
with our previous result, eq. ^. Table ^ shows that all three fits lead to stable values of rrf"^ 
as well as k^^^ and Kgf^*. This gives us confidence in our method of choice: we employ a 
semi-constrained fit to extract the strange quark mass. 

The data and used as follows: 



Determine k''*i'^'^s° from rrisv by matching 



1.16 , (18) 



where i^^i^^ is given by eq. ^. Alternatively, the K can be used^ 
• Determine i^^'^'^^'^S'^ from ttt-sv by matching 

TOi/ fK^'S'^* ^strange N 

mvuMf^sc. j = ^^ = 1.326 . (19) 

/ light \ ]\/r ^ 

my2_ss('«sea ) ^^^P 

Table ^ shows the results for all three particles employed. Good agreement is found between 
the (j) aiid the K*, whereas the K favours a slightly larger ft;Sti'ange_ 'pj-^jg jg illustrated in figure 
0, where ^strlngc is plotted as a funtion of 

To check for a systematic error we have also tried constrained fits to the vector particle. 
Encouragingly, we find only negligible change (see table ^). 

Our final value is: 

^strange ^ 0.15615(20)^*'* (20)^^^* , (20) 
where the systematic error covers the spread for all three particles. The conversion yields: 

m^''^'(2 GeV) = 140(20) MeV . (21) 

Before we turn to a discussion of our results, we comment on a method which has been proposed 
for the extraction of /^strange ^^^i^ 2 dynamical quarks (see 0, for example). It circumvents 
the introduction of valence quarks unequal to the dynamical quarks and determines n^^^^^s^ 

®In which case we match 
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6.6 



6.5 



6.4 



6.3 



6.2 




6.30 



6.35 



6.40 



Figure 3: Chiral behaviour of ^stringc ; as determined from K*{0), 0(0) and the K{0). The 
strange quark mass is calculated with the light quark mass fixed to its physical value. At the 
very left we plot the value for a strange quark in a strange sea. Lines are calculated from eq. ^ 
and eq. ^ 



particle 


^strange constrained 


^strange gemi-constrained 


K 


0.1564 + I 


0.1565 t 2 


K* 


0.1561 t 2 


0.1559:1: 2 


4> 


0.1561 t 2 


0.1559:^ J 



Table 3: Collection of results for K^trange^ 



by matching the ratio 

8 



..strange 



J ^ using the symmetric fit only. Adopting this procedure, one finds 

0.15709 Zil and m^''^''(2 GeV) = 80(8) MeV; this point is shown on the very left in 
figure ^. However, this method implies that the (j) mass is measured with strange valence quarks 
in a sea of strange quarks. One would prefer the mass to be determined in a sea of light quarks, 
at which, naively, one expects the effects of dynamical sea-quarks to be more pronounced. It is 
interesting to see that the effect of light sea-quarks is fairly significant: the 0, when measured 
in a sea of strange quarks, yields the largest of all estimates for ft;Strange_ 



4 Discussion 

Our results 



r4g(2GeV): 

strange /n /- \/\ 



2.7(2) MeV 
: 140(20) MeV 



2 dynamical light quarks 



(22) 
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1/K 



1/K = 1/K^ 

sea V 



'1/K, 



1/K. 



1/K„ 



1/K. 



1/K 



l,v 



Figure 4: The values of -K- and —h^ are plotted versus The situation is displayed schemat- 



ically around the region of the symmetric line k^, 



are to be compared to the analogous quenched values at corresponding /^quenched = 6.0: 



r^f (2GeV): 

^sgnge(2GeV) 



5.5(5) MeV 
: 166(15) MeV 



quenched 



(23) 



Errors due to the finite lattice spacing and the finite volume are not included in (22) and (|2^). 

Compared to the quenched result, which is in good agreement with previous lattice calcula- 
tions we observe a much smaller dynamical light quark mass, whereas the strange masses are 
compatible within errors. Our dynamical result for the quark mass ratio, TT^strange^j^iight ^ 
differs significantly from chiral perturbation theory estimates|^] or sum rule resultsQ. This may 
be due to the fact that the effects of strange sea quarks can only be partially accounted for in 
an Nf = 2 simulation. 

We note that an extrapolation of quenched quark masses to a = yields much smaller values 
than eq. |2|, m"g'^* = 3.4(5) and m'^*''^'^^^ = 100(23) @]. Nevertheless, the ratio rn''^^'^''^'^ /m^is^^ 
from quenched simulations is fairly independent of the lattice spacing a. This may be different 
in full QCD, as cutoff effects could show up differently in sea and valence quarks. Dynamical 
results at other couplings are needed before one can compare results in the a = limit. 

To understand the drop in the mass of the light quark from A'^^ = to A'^^ = 2 at fixed 
coupling we have analysed our dynamical data in a manner suggested by quenched lattice sim- 
ulations. To this end, we have defined a quark mass at fixed sea-quark through: 



m 



-(—--) 



(24) 



Setting m|« ^^(k! 



clt f^'\J 



7^ Kggg^ effectively forces an unphysical pion to become massless. 



Using eq. 16 we can calculate the quenched critical kappa at fixed sea-quark as: 



ci + 



C3 



C2 



(25) 



'^v('^sea) 

to be compared to the true critical kappa ^j— = — ^^^^^ ■ By construction (eq. the values 
of — r will lie on a straight line hitting on the symmetric line. The condition = -4— 
has two solutions, C3 = and 03 + 02 = — ciKsea- The latter case is identical to the true critical 
kappa. The trivial condition C3 = corresponds to vanishing sea-quark dependence. 
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Proceeding as in quenched simulations and measuring a bare light quark mass at each of the 
three sea-quark values we find light quark masses very similar to that of the quenched simulation 
(5.7(4), 5.6(3), 5.4(3) MeV) with only a slight downward trend. An extrapolation of these quark 
masses in the sea-quark to the critical point would yield the value A2 in figure |^, while the true 
value is Ai; the latter represents the value of quark masses for a physical pion in a physical 
sea-quark, whereas A2 is given in a sea of massless quarks. However, if we try to repair A2 by 
extrapolating to the light sea-quarks we have to give up working at the physical pion mass since 
then the critical kappa ^ is too low. This means that the light quark mass, which we wish 
to obtain from the physical pion mass and in a sea of physical up and down quarks, cannot be 
obtained by extrapolating values obtained at fixed sea-quark to either the critical or the light 
sea-quark mass. Figure Q also illustrates that guenc/ied Wilson-like light quark masses away from 
should be treated with caution since they lead to negative quark masses when measured 

sea 

with respect to the physical point -j— , (-4 — -J—) < 0. 



Finally, we note that a study of finite size effects is under way|ll]. We shall include an 
additional sea-quark value which will allow us to study the effects of higher order terms in chiral 
perturbation theory. For the simulation presented here, we postpone such a discussion to 0. 
An analysis of the bottomonium spectrum, currently in progress, will allow us to use a lattice 
spacing obtained from the 15 — IP splitting, which should be less sensitive to lattice artefacts. 
Much more study is needed, of course, to gain control over these. 
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